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Multi-orbital Hubbard model Hamiltonians for the undoped parent compounds of the Fe-pnictide 
superconductors are here investigated using mean- field techniques. For a realistic four-orbital model, 
our results show the existence of an intermediate Hubbard U coupling regime where the mean-field 
ground state has spin stripe magnetic order, as in neutron scattering experiments, while remaining 
metallic, due to the phenomenon of band overlaps. The angle-resolved photoemission intensity and 
Fermi surface of this magnetic and metallic state are discussed. Other models are also investigated, 
including a two orbital model where not only the mean- field technique can be used, but also Exact 
Diagonalization in small clusters and the Variational Cluster Approximation in the bulk. The 
combined results of the three techniques point toward the existence of an intermediate-coupling 
magnetic and metallic state in the two-orbital model, similar to the intermediate coupling mean- 
field state of the four-orbital model. We conclude that the state discussed here is compatible with 
the experimentally known properties of the undoped Fe-pnictides. 

PACS numbers: 71.10.Fd, 72.10.Di, 72.80.Ga, 79.60.-i 



I. INTRODUCTION 



The discovery of superconductivity in the Fe pnictides 
has opened an area of research that is attracting consid- 
erable attentionii^i^i^i^i^iii^ In the early stages of these 
investigations, the layered structure of the Fe pnictides 
superconductors,— i^i^i^i^i^i^i^ the existence of a magnetic 
spin striped state revealed by neutron scattering in the 
undoped limit,— and their large superconducting crit- 
ical temper at ure^i^i^i^i^i^i^i^ motivated discussions on a 
possible close relation between these new Fe-based mate- 
rials and the high-temperature cuprate superconductors. 
However, it was clear from the initial investigations that 
there were substantial differences as well: for example, 
the resistivity vs. temperature curves of the parent com- 
pounds^ ^^i^i^i^i^i^ii^ do not show the characteristic Mott 
gapped behavior of, e.g., LaCu04. In fact LaOFeAs be- 
haves as a bad metal or semiconductor ^i^i^i^i^i^i^i^ but 
not as an insulator. Moreover, the magnetic moment 
in the spin striped state of LaOFeAs is much smaller 
than expected i^i^ Although further neutron scattering 
research has shown that the magnetic order parame- 
ters are larger in other Fe pnictides, their values are 
still below those anticipated from band structure calcu- 
lationsi^i^ii^i^ii^ or from the large Hubbard U limit of 
model Hamiltonians^^ (unless couplings are in a spin frus- 
trated regime^). In summary, the parent compounds 
of the Fe superconductors behave in a manner different 
from the parent compounds of the Cu-oxide supercon- 
ductors because the zero temperature resistivity is finite 
and the magnetic order weak. However, the pnictides 
are also different from BCS materials, where the normal 



state is a non- magnetic metal with low resistivity. Then, 
the Fe superconductors appear to be in an intermediate 
regime of couplings, somewhere in between, e.g., MgB2 
and the Cu oxide superconductors ^^^^ The "antiferro- 
magnetic metallic" nature of LaOFeAs is clearly different 
from an antiferromagnetic insulator or a non-magnetic 
metal. Further confirming the need to focus on the in- 
termediate coupling regime for the Fe pnictides parent 
compounds, a pseudogap has been observed— i^i^i^ in 
their density of states (DOS), which is different from the 
featureless DOS of a good metal or the gapped DOS of 
an insulator. 

Much of the current theoretical ef- 

f Q^^^lS ^25 ^26 ^27 ^28 ^2 9 ^30^ 31 ^ 32 ^ 33 ^34 ^35 ^36 ^^38 ^^40 ^41 ^42 ^43 ^44 ^45 ^46 

has focused thus far on two well defined limits. On one 
hand, band structure calculations have reported the 
existence of Fermi pockets,— i^i^i^i^ which were con- 
firmed by photoemission experiments i^^i^^i^^i^^i^^i^^i^^ 
On the other hand, intuition on the physics of the model 
Hamiltonians is often gained by investigating the large U 
regime. Although results by some of us for a two-orbital 
model using numerical techniques have already shown 
that a small magnetic order can be accommodated 
at intermediate couplings^ this range of couplings is 
typically the most difficult to handle using computer 
simulations and, moreover, more bands are expected to 
be of relevance for a better quantitative description of 
the Fe pnictides. 

In this publication, multi-band Hubbard models are 
investigated using mean- field and numerical techniques. 
Our most important result is the discovery of an inter- 
mediate Hubbard U coupling regime where the ground 
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state is an antiferromagnetic metal. More specifically, 
for U larger than a critical value Ud the spin striped or- 
der develops with continuity from zero. In spite of a gap 
at particular momenta, the overall state remains metallic 
(there is a nonzero weight at the chemical potential in the 
DOS) due to the phenomenon of band overlaps. Further 
increasing the coupling to Uc2^ a fully gapped insulator 
is stabilized. Thus, the intermediate regime Uci<U<Uc2 
is simultaneously (i) magnetic with a small order param- 
eter and (ii) metallic. Moreover, the DOS reveals the 
existence of a pseudogap in this regime. All these prop- 
erties are compatible with our current knowledge of the 
Fe pnictides parent compounds. 

While these results are interesting, they were obtained 
using mean- field approximations. For the realistic case 
of four or five orbit als, it is difficult to obtain reliable 
numerical results to confirm the mean- field predictions. 
However, for the case of two orbit als, calculations can 
be carried out using simultaneously the mean-field tech- 
nique, that also reveals an intermediate coupling regime 
similar to that of the four and five orbit als models, to- 
gether with the Exact Diagonalization (ED)^ and Vari- 
ational Cluster Approximation (VCA) ^^i^^ methods used 
before. Below, it is reported that the results using these 
computational methods are compatible with those of the 
mean-field approximation, providing confidence that the 
mean-field method may have captured the essence of the 
problem. However, it is certainly desirable that future 
investigations confirm our main results. 

For completeness, here some related previous litera- 
ture is briefly mentioned. The band overlap mecha- 
nism for an insulator to metal transition has been exten- 
sively studied before using band structure calculations 
in a variety of contexts, such as TlCl and TlBr,^^ solid 
hydrogen^i^ and bromine under high pressure. Closer 
to the present results, the existence of an intermedi- 
ate Hubbard U regime with an antiferromagnetic metal- 
lic state was previously discussed by Duffy and one of 
the authors (A.M.)^^ in the context of the high tem- 
perature Cu-oxide superconductors and using the one- 
band Hubbard model, after introducing hopping terms 
f and between next-nearest-neighbor sites. Den- 
sity functional methods were also used before to dis- 
cuss AF-metallic states and band-overlap insulator-metal 
transitions i^i^ Within dynamical mean- field theory, an 
AF-metallic state has also been discussed. Experimen- 
tally, itinerant antiferromagnetic states were found in the 
pyrite NiS2-xSeajr^i^ in heavily doped manganites,^^ in 
ruthenates,^^ in organic conductors'*^ and in several 
other materials. An incommensurate spin density wave 
was also reported in metallic y2-y03^ The results dis- 
cussed in this manuscript establish an interesting con- 
nection between the Fe pnictides and the materials men- 
tioned in this paragraph. 

The organization of the paper is as follows. In Section 
II, results for a four-orbital model are presented. This 
includes a discussion of the model, the mean-field tech- 
nique, and the results, with emphasis on the intermediate 



coupling state. The photoemission predictions for this 
state are discussed. In Section HI, we show results for the 
two-orbital model using ED^ and VCA method^*^ in 
addition to mean- field approximations. Section IV con- 
tains our main conclusions. 



II. RESULTS FOR A FOUR-ORBITAL MODEL 

In this section, we will describe a possible minimal 
four-orbital model for the Fe-based superconductors. 
This model presents a Fermi surface similar to that ob- 
tained with band structure calculations. 



A. The four-orbital model 

Previous studies have suggested that the Fe-As planes 
are the most important substructures of the full crystal 
that must be analyzed in order to reproduce the physi- 
cal properties of the Fe pnictides close to the Fermi sur- 
face. We consequently focus on these Fe-As planes in the 
present study. The effective Fe-Fe hopping Hamiltonian, 
using As as a bridge, can be obtained within the frame- 
work of the Slater-Koster (SK) formalism. By this pro- 
cedure, here we will construct a minimal model defined 
on a Fe square lattice, consisting of the four Fe d or- 
bitals xz, yz^ xy^ and — y^- It is assumed that the 
dsz'^-r'^ orbital lies at a substantially lower energy and is 
thus always filled with two electrons ^^i^i^i^i^ While 
the SK procedure is not as quantitatively accurate as 
a full band-structure calculation, it can still provide the 
proper model Hamiltonian, because it correctly takes into 
account the geometry of the system and illustrates which 
orbitals are connected to one another at different lattice 
sites. Thus, our procedure here will be to use the SK 
method to construct the formal model, and then obtain 
the actual numerical values of the hopping parameters 
via comparison with band-structure calculations. 

The Hamiltonian H = Hq + i^int includes two parts: 
the hopping term Hq and the interaction term i^int • The 
hopping term in real space reads 

^0 = E E E (n^"4.,.^j-.- + f^-^) ' (1) 

i,j cr 

where dl ^ ^ creates an electron at site i with spin a 
on the /i-th orbital (/i = 1,2,3,4 stands for the xz, 
yz^ xy^ and x'^-y'^ orbitals, respectively). Here, hop- 
pings at nearest-neighbors (NN) and also at next-nearest- 
neighbors (NNN) along the plaquette diagonals were con- 
sidered. The hopping tensor T^^j^ has a complicated real- 
space structure that will not be reproduced here. Hq has 
a simpler form when transformed to momentum space: 

^0 = EEE^'''''(k)<M,.'^k,.,., (2) 
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with 

T^^ = - 2^2 cos /c^ - 2^1 cos /c^ -4^3 cos /c^ cos A;^, (3) 
T^^ = -2ticos/c^ -2t2Cos/c(^ -4t3Cos/c^cos/c' (4) 



= -4^4 sin A;; sin (5) 
T^^ = -2t5(cos(A;; +7r) +cos(A;(^ +7r)) 

-4^6 cos(A;; + tt) cos{ky + tt) + A^y, (6) 

T^^ = -4zt7sin/c; + 8it8sinA;;cosA;(^, (7) 

T^^ = -4it7 sin ky + Sits sin ky cos fe^ , (8) 
T"^^ = -2ti7(cos(/c; +7r) +cos(/c(^ +7r)) 

-4^9 cos(/c^ + tt) cos{ky + tt) + A2,2_^2 , (9) 

_ -4ztiosinA;^, (10) 

T^^ ^ 4ztiosinA:;, (11) 

T^^ = 0. (12) 



Equation ([2]) represents a matrix in the basis {dl. ^ ^} 
where k = kMf /i = 1 or 2 and k = + Q if /i = 3 or 
4 with — TT < kx^ky < TT, k' is defined in the reduced 
Brihouin zone corresponding to the two-Fe unit ceh, and 
Q = (tt, tt). In other words, the above expressions couple 
states with momentum k' for orbitals xz and yz to states 
with momentum k' + Q for orbitals xy and x^ — y^ . The 
momentum Q appears after considering the staggered lo- 
cation of the As atoms above and below the plane defined 
by the Fe atoms. As already mentioned, the mathemat- 
ical form of this model arises directly from the Slater- 
Koster considerations. This problem is equivalent to an 
eight-orbital model with a Hamiltonian expanded in the 
basis {d\^, ^ ^, d\^,j^Q^ /x ai- As a result, the eight-orbital 
band structure and Fermi surface (Fig. [2]) are obtained by 
"folding" the results in the equivalent four-orbital prob- 
lem (Fig.[l]). 

The actual values of the hopping parameters could in 
principle be obtained from the overlap integrals in the 
SK formalism. ^^'^^ However, to properly reproduce the 
Fermi surface obtained in the Local Density Approxima- 
tion (LDA)i^J^iiiii^J£ it is better to fit the values of 
those hoppings. The parameters used, as well as the on- 
site energies for the xy and x^ — y^ orbitals, are listed 
in TableU The on-site energy term is given by ^ A^n(^ 
(standard notation) and it is part of the tight-binding 
Hamiltonian. 

In Figs. [Hand [21 we show the band structure and the 
corresponding Fermi surface for the four-orbital tight- 
binding Hamiltonian in Eq. ([2|), using Table [H Since 
we assume that the Zz^ — orbital is always doubly 
occupied, the chemical potential in the undoped case is 
determined by locating n = 4 electrons per site in the 
four bands considered here. As shown in Fig. [ijb), two 
hole pockets centered at (0,0) (arising from the ol\ and 
Oi2 bands) and four pieces of two electron pockets cen- 
tered at (0, tt) and (tt, 0) (from the f5\ band) are obtained. 
The shape of the Fermi surface qualitatively reproduces 
the band-structure LDA calculations^i^ii^ii^ii^ after a 



TABLE I: Fitted hopping parameters and on-site energies for 
the four-orbital model used in this section (in eV units). 
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FIG. 1: (Color online) (a) Band structure corresponding to 
the four-orbital tight-binding Hamiltonian Eq. (|2|, using the 
LDA fitted values of the hopping parameters provided in Ta- 
ble [T] The Fermi surface is formed by two hole-like bands (ai, 
0^2 ) and one electron- like band (/3i). The chemical potential 
is at 0. (b) The topology of the corresponding Fermi surface. 

45^ rotation about the center of the First Brillouin Zone 
(FBZ), as presented in Fig.[2jb), due to the rotation from 
the Fe-Fe axis to the Fe-As axis. 

To study the relation between the orbital hybridization 
and the Fermi surface topology, the projected weight of 
each orbital at both the hole and electron pockets were 
calculated. These weights are defined via the eigenvec- 
tors of i^o: ^/i,A(k) = lEa l^k,/i,(7;Ap, whcrc A de- 
notes the band index (ai, 0^2, A, /^2), and /i refers to 
the four d orbitals. The matrix t/k,/i,cr;A diagonalizes 
the system (see Eq. ([T8|) below). An example of the 
angle-resolved weights in momentum space are shown 
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FIG. 2: (Color online) (a) Band structure of the eight-orbital 
problem in the reduced BZ obtained by "folding" the results 
presented in Fig. [TJa). (b) Fermi surface of the eight-orbital 
problem obtained by folding the FS obtained in Fig. [2b). 
The FS in the first (second) BZ are indicated by continuous 
(dashed) Unes. 

in Fig. [3l The two hole pockets centered at (0, 0) 
mostly arise from the xz and yz orbitals, compatible 
with LDAi^ii^ii^ii^ii^ and with much simpler descrip- 
tions based only on two orbitals J^i^ The electron pocket 
centered at (tt, 0) ((0,7r)) arises mainly from the hy- 
bridization of the xz {yz) and xy orbitals (not shown). 
These results are also qualitatively consistent with those 
from the first-principles calculations However, there 
are some quantitative discrepancies that lead us to be- 
lieve that probably longer-range than NNN plaquette- 
diagonal hoppings are needed to fully reproduce the LDA 
results including orbital weights. Nevertheless, the dis- 
cussion below on the metallic magnetic phase at interme- 
diate couplings is robust, and we believe it will survive 
when more complex multi-orbital models are used in the 
future. 

Note that the eigenenergies (band dispersion) along the 
(0, 0) (tt, 0) and (0, 0) (0, tt) directions are symmet- 
ric about (0, 0), but the eigenvectors (W^^x) show a large 
anisotropy. For instance, at the Fermi level the ai band 
is almost xz-like along the (0, 0) (tt, 0) direction but 
almost yz-\ike along the (0, 0) (0, tt) direction. Be- 
low, it will be discussed how this anisotropy affects the 
mean-field results for the interacting system. 




0/71 

FIG. 3: (Color online) The projected orbital weight W^^x of 
states at the Fermi surface. Shown, as example, are results 
for the outer hole pocket centered at (0, 0). The definition of 
G is given in the inset. 

Let us now consider the interaction term,^^ which reads 
Hint = ^^ni,^,|ni,^,| + ([/'- ^) ^ ni,^ni,^ 

-2J ^ Si,^ • Si,,, (13) 

where Si,^ is the spin (charge density) of orbital fi 

at site i, and ni^^ = ni^^^|+ni^^^|. The first term is a Hub- 
bard repulsion for the electrons in the same orbital. The 
second term describes an on-site inter-orbital repulsion, 
where the standard relation U' = U — J/2 caused by rota- 
tional invariance is used.^^ The last term in Eq. (p!3|) is a 
Hund term with a ferromagnetic coupling J. A complete 
description would also require a pair-hopping interaction 
similar to the last term of Eq. ([24j) . where the interac- 
tion term for the two-orbital model is shown. But ED 
was used to test its impact in the case of two orbitals, 
and it was not found to be important. Consequently, it 
was neglected in the mean field treatment. 

B. The mean-field approach 

To study the ground state properties of the system, we 
apply a mean-field approximation to the model Hamil- 
tonian described by Eqs. ^ to (fT3|) . We follow here 
the simple standard assumption of considering only the 
mean- field values for the diagonal operators )2i 

where q is the ordering vector of the possible magnetic or- 
der, and are mean-field parameters describing the 
charge density and magnetization of the orbital /i, and 
the rest of the notation is standard. Applying Eq. (p!4|) 
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to i^int, the mean- field Hamiltonian in momentum space 
can be written as 

HmF = Hq -\- C -\- ^ ^/i<^k,/i,cr^k,/i,cr 
k,/2,o" 

+ Yl '?P,<T(4,^,a^^k+q,p,a +4+q,^,a^^k,^,<T), (15) 
k,/i,cr 

wiiere k runs over tiie extended FBZ, Hq is tiie iiopping 
term in Eq. ([2]), 

C= - NU^^fnl- ^ruA - N{2U' - J)Y, n,n, 
NJ 

is a constant, N the lattice size, and we used the defini- 
tions 

= Un^ + (2C/' " E (16) 

Vfi,(T = h7m^ + J ^ . (17) 

The above mean-field Hamiltonian can be numerically 
solved for a fixed set of mean-field parameters using stan- 
dard library subroutines. The parameters and are 
obtained in a self-consistent manner by minimizing the 
energy. In practice an initial guess for and serves 
as a set of input parameters for a given value of the cou- 
plings U and J. The mean-field Hamiltonian is then diag- 
onalized and and are reevaluated using Eq. (p!4|) . 
This procedure is iterated until both and have 
converged. During the iterative procedure Ti^=A was 
enforced at each step, such that the total charge density 
is a constant. Note that in the mean-field approximation 
an electron with momentum k is coupled to an electron 
with momentum k + q, where q is the vector associated 
with the magnetic ordering. Then, the Hamiltonian in 
Eq. (p!5|) is solved only in the magnetic reduced Brillouin 
Zone with only half of the size of the unfolded FBZ. 

The numerical solution of the mean-field Hamiltonian 
immediately allows for the evaluation of the band struc- 
ture, the density of states (DOS), and the magnetization 
(M = m^) at the ordering wavevector q. Moreover, 
we can also calculate the photoemission spectral func- 
tion. Assuming that the mean-field Hamiltonian Eq. (p!5|) 
is diagonalized by the unitary transformation 

<^k,/i,a = ^ ^k,/i,(7;A7A, (18) 
A 

i^MF = ^PatItA, (19) 
A 

then the spectral function is given by 

A(k, '^) = Y.Y. t^k,p,.;Af/l;k,.,.^(^ - pa). (20) 

A M,Z^,cr 



In practice ^(cj— pa) was here substituted by ^ ^s^^JL^^p 
with a broadening e = 0.05 eV. 

C. Mean- field results 

In this subsection, the mean-field results for the four- 
orbital model previously described are presented. Exper- 
imentally, a (bad) metallic phase with spin-stripe mag- 
netic order at wavevector s (0, tt) was observed in the un- 
doped compound LaOFeAs.^^^^ It is then important to 
investigate in the mean-field approximation the proper- 
ties of such a spin striped ordered state. Here we show 
the numerical results for the mean-field approximation 
defined on a 100 x 100 square lattice. 



1. Magnetic order 

To study in more detail the spin striped ordered state, 
the ordering wavevector q is assumed to be (0,7r) in the 
mean-field approximation. In Fig. HJa), the evolution of 
the magnetization vs. /7, at J=U /4, is shown. At a crit- 
ical value Uci = 1.90 eV, the (0, tt) order— starts to grow 
continuously from zero. This "stripe" magnetization in- 
creases slowly until it reaches Uc2 = 3.75 eV, where it 
changes discontinuously, showing the characteristic of a 
first-order transition (see discussion for the origin of this 
transition later in this section). Note that for U ^ Uc2^ 
the stripe magnetization for the (0, tt) state is smaller 
than 1.0 (with a normalization such that the maximum 
possible value is 4.0), indicating that there is an interme- 
diate U regime that can accommodate the rather weak 
striped magnetic order found in the neutron scattering 
experiments. ^^^^ In Fig. [41(b), magnetization curves for 
the same state at various values of J are presented. Two 
transitions for all J /U ratios studied are observed, simi- 
larly as for J = [7/4 in (a). 

2. Band structure and Fermi surfaces 

Let us analyze in more detail the spin striped ordered 
state. Since varying J/U does not significantly alter the 
results, the value J = U/A\s adopted in the rest of the 
analysis. In Figs. O [6l and [3 the band structures and 
Fermi surfaces of this state at several values of U are 
extracted from the calculated mean-field photoemission 
spectral function data. Both the band structure and the 
Fermi surface at the critical point Ud = 1.90 eV are, of 
course, identical to those at /7 = in Fig. [H This non- 
interacting electronic description of the system changes 
gradually upon the establishment of the striped magnetic 
order. As discussed in more detail below, gaps open 
at particular momenta, while other bands crossing the 
Fermi surface do not open a gap. Thus, this is a metallic 
regime with magnetic order. Finally, at the magnetiza- 
tion discontinuity a full gap develops. 
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FIG. 4: (Color online) (a) The mean-field evolution with U of 
the magnetization of the state ordered at wavevector (0,7r), 
with J — U/4:. Shown are the two critical values: one where 
the magnetization becomes nonzero and a second one where 
a discontinuous behavior is concomitant with a metal to in- 
sulator transition, as described in the text, (b) The same 
magnetization curve as in (a), but for several values of J/U. 
The magnetization shown in (a) and (b) is normalized such 
that the maximum value is 4, corresponding to four spin po- 
larized electrons, one per orbital, at each site. 



3. Bands of magnetic origin in the spin striped state 

Due to the off-diagonal term in Eq. (fT5|) . an electron 
with momentum k is coupled to another with momentum 
k + q if their orbital characters are the same. This gener- 
ally leads to avoided level crossings and opens a gap pro- 
portional to the magnetization. If the hole and electron 
Fermi surfaces can be connected by the vector q, then a 
magnetic state with spin ordering at q is stabilized over 
the non-magnetic state due to the Fermi surface nesting 
effect. The magnetic state gains energy via the open- 
ing of a gap near the Fermi level. This magnetically or- 
dered state significantly changes the band structure of the 
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FIG. 5: (Color online) Mean-field photoemission band struc- 
ture of the spin- striped (0,7r) state in the energy window [- 
2 eV, 2 eV] at [/ = 1.90 eV, V = 2.50 eV, V = 3.20 eV, 
and U = 4.00 eV, from top to bottom, with J^U/A. The 
first U is at the first transition, and still has the shape of the 
noninteracting limit. The second two values of U are in the 
intermediate coupling regime, and the presence of states at en- 
ergy (location of the chemical potential) indicate a metallic 
state. The last coupling, U = 4.00 eV, is above the second 
critical point, and thus already in the gapped regime. In the 
two intermediate couplings, 2.50 and 3.20 eV, weak bands not 
present in the non- interacting limit (top panel) are revealed. 
These are the bands caused by the striped magnetic order, 
which should be observable in photoemission experiments. 



non-interacting case, opening gaps and giving rise to the 
emergence of new bands of magnetic origin, sometimes 
called the "shadow bands", as shown in Fig.[5](b) and (c), 
and with more detail for a special case in Fig. [6l Similar 
issues were discussed in the context of high temperature 
superconductors, where the existence of bands caused 
by the staggered magnetic order were extensively stud- 
ied before In fact, experiments for undoped cuprates 
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FIG. 6: (Color online) Mean-field photoemission band struc- 
ture of the spin- striped (0,7r) state in the energy window [- 
0.8 eV, 0.8 eV] at [/ = 2.50 eV. Bands caused by magnetic 
order that are not present in the non-magnetic case U = 0.0 
are here shown in more detail than in Fig. [5] 



revealed a photoemission spectral function in excellent 
agreement with theoretical expectations^ i.e. contain- 
ing the predicted bands generated by the magnetic order. 
Thus, it is to be expected that the spin striped order 
of the Fe pnictides should also produce magnetically- 
induced bands in the undoped limit, and even in the 
doped case if the magnetic correlation length remains 
large enough. 



4. The first order transition at Uc2 

Regarding the second transition at /7c2, the qualita- 
tive reason for its presence lies in the incompatibility of 
the intermediate coupling metallic state with the large U 
limit. The mechanism of nesting that causes the special 
features of the intermediate U state previously discussed, 
including the survival of portions of the Fermi surface, 
will lead to an energy that eventually cannot compete 
with a fully gapped state at large U at the electronic 
density considered here, thus a transition must eventu- 
ally occur. 

But why is the second transition discontinuous? In 
Fig.[5]the coupling between the hole pockets and the elec- 
tron pocket at (0, tt) leads to the distortion of the Fermi 
surface in the spin-striped state. Such a coupling between 
states with a specific orbital symmetry and momentum 
also accounts for the metallic nature of the spin-striped 
state: the electron pocket at (tt, 0) is almost undistorted 
for U < Uc2' However, when U is approaching the sec- 
ond critical value Uc2^ the peak at (7r,7r) with occupied 
states in the ai band becomes energetically closer to the 
Fermi level (its energy is increasing with U). If the Fermi 
surface nesting effect could be neglected at /7c2, then a 
smooth behavior would be observed since the charges 
could transfer continuously from the peak at (tt, tt), after 
crossing the Fermi level, to the peak at (0,0). However, 
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FIG. 7: (Color online) Mean-field photoemission Fermi sur- 
faces in the spin-striped (0,7r) state at [/ = 1.90 eV, U = 
2.50 eV, and U = 3.20 eV, from top to bottom. The re- 
sults are obtained via A(k, 00) using a window of 20 meV cen- 
tered at the Fermi energy. Shown are results obtained from 
an equal weight average of data using A(k,a;) and A(p,a;), 
where \i={kx,ky) and p={kx, —ky). By this procedure the 
results are properly symmetrized under rotations in the non- 
magnetic phase. The anisotropy of the results in the spin- 
striped intermediate U region do appear because the spin or- 
der breaks rotational invariance. 



note that the valley of Pi band at (tt, 0) and the peak of 
the ai band at (tt, tt) have both a partial xy symmetry. 
Moreover, they are connected by the vector q = (0,7r). 
Thus, it will be expected that a gap close to the Fermi 
level will open to minimize the energy. At Uc2^ the sys- 
tem gains maximal energy by opening a finite gap at both 
(tt, 0) and (tt, tt), and lowering the energy of the a bands 
at (0, 0) such that they become fully occupied. This leads 
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to discontinuous changes in the population of the indi- 
vidual orbit als, producing discontinuous changes in the 
orbital magnetizations and, concomitantly, a finite gap. 
Since in the real undoped Fe-pnictide materials the full- 
gap regime is not realistic, then we can proceed with the 
rest of the analysis below without further consideration 
of this discontinuity in the magnetization. 



5. Anisotropic Fermi surface 
in the undoped parent compound 

The appearance of spin-striped magnetism with 
magnetically-induced bands leads to an anisotropic dis- 
tortion of the Fermi surface in the magnetically ordered 
state. For instance, consider the q=(0,7r) mean-field 
state. To predict the pattern of gaps in this state, the 
information previously discussed in, e.g.. Fig. [3] is im- 
portant. Through the projected weights, it was observed 
that the j3i electron pocket centered at (0, tt) has mainly 
yz-\-xy symmetry, while the (tt, 0) electron pocket is 
mainly xz-\-xy. A dominance of the yz symmetry is found 
in the ai hole pocket along the (0, 0) (0, ±7r) direc- 
tions, and in the 0^2 hole pocket along the (0, 0) (±7r, 0) 
directions. Then, for the magnetic state with q=(0,7r), 
a gap should open in the (0, 0) (0, ±7r) directions 
for the ai (inner hole pocket) band and for the elec- 
tron pocket at (0,7r), but the 0^2 (outer hole pocket) 
band remains gapless. In addition, a gap opens in the 
(0, 0) (iTT, 0) directions for the a2 band, but both the 
inner hole pocket and the (tt, 0) electron pocket remains 
gapless. These results are indeed observed numerically 
as shown in Figs. [5l^b,c) 

Due to these anisotropic gaps, the Fermi surfaces of the 
a bands change their topology from two pockets Fig.fTJa) 
to four arcs that are very close to one another, as shown 
in Fig. [TJb). In between the arcs, there is actually a 
nonzero but weak intensity at the location of the original 
hole pockets. For LaOFeAs, since a weak striped order 
has been observed, the Fermi surface is expected to have 
a similar shape as that shown in Fig. [T^b). However, 
topology of the Fermi surface in this state is sensitive to 
the value of the coupling U : further increasing this cou- 
pling the four arcs merge into a single pocket Fig. EJc). 
Then, ARPES experiments can provide valuable infor- 
mation about the exact shape of the Fermi surface and, 
thus, the interaction strength. 

Besides the anisotropic gaps near the Fermi level, gaps 
far from the Fermi level also exists (see, for instance. 
Fig. [51(b)). This cannot occur in a one-band model be- 
cause the opening of these gaps would not provide any 
energy gain. However, such gaps are possible in more 
complex multi-orbital models: since the off-diagonal term 
in Eq. (p!5|) depends on the magnetization that is con- 
tributed by each orbital, once the magnetic state is sta- 
bilized by the opening of a gap near the Fermi level, the 
off-diagonal term becomes non-zero even for the bands 
far from the Fermi level. 




FIG. 8: (Color online) The angular resolved weight W' at 
[7=0 obtained using Eq. (|2ip for the two a hole-like pockets. 
The definition of G is the same as in Fig. [3l 



A rather surprising result is that the intensities of 
the spectral function of the two a pockets display an 
anisotropy even in the non-magnetic phase (not shown). 
This is puzzling since from Fig. [3l it can be observed 
that Xl^VK/2,A(k) = 1 for any k. Naively, this would 
lead to an isotropic A(k^Lj). However, from Eq. ([2Q|) . 
A{k,u,) = - px)W;,{k), where 



W',{k) = J2\J2u^,^,..,x\'. (21) 



Here, note that an anisotropy could arise from the inter- 
ference between different orbit als. To observe this more 
explicitly, W'^ for the two a hole pockets at /7 = is 
shown in Fig. [8l These functions are not constant but 
oscillate between and 2. From Fig. [HI it is observed that 
for the ai band W reaches its maximum at, e.g., — 7r/4, 
corresponding to the (0,0) (— 7r,7r) direction, while 
for the a2 band, the maximum of W is at, e.g., 7r/4, i.e., 
in the (0, 0) (tt, tt) direction. Such an anisotropy in 
74(k, uo) is not a consequence of the mean-field approxima- 
tion, but is related to the multi-band nature of the model 
itself. Interestingly, such an anisotropy in the topology 
of the Fermi surface in the A(k, uo) data may account for 
the anisotropic features of the hole pocket in the recent 
ARPES experiment on LaOFeP^ where a long-ranged 
striped magnetic order is absent.— However, in Fig. [7] 
this problem is avoided by a symmetrization procedure 
(see caption of Fig. [7j) that restores rotational invariance 
to the non-magnetic state. Thus, the anisotropy of the 
important intermediate U regime fully originates in the 
lack of rotational invariance of the spin striped state with 
a wavevector (0,7r) or (7r,0). 
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6. Metallic magnetically ordered phase at intermediate 
couplings and existence of a pseudogap 

As already remarked, Fig. [5] indicates that for moder- 
ate U the magneticahy ordered system is stih gapless, i.e. 
it is in a metahic phase with a finite Fermi surface due 
to the phenomenon of "band overlapping" described in 
the Introduction. This intermediate state is also revealed 
via the evolution of the DOS in Fig. [9l As displayed in 
Fig. [Ha), a pseudogap in the DOS near the chemical 
potential exists in the regime Ud < U < Uc2- If the 
availability of states at the Fermi level is assumed to be 
directly related with transport properties, a DOS pseu- 
dogap suggests bad-metallic characteristics in the inter- 
mediate U regime. The pseudogap turns into a hard gap 
at /7 > /7c2, where the system becomes an insulator. In 
Fig. [H^b) the value of the DOS at the chemical potential 
vs. U is plotted. Two transitions can be easily identi- 
fied. For U < Uci^ A^(m) is a constant. It decreases con- 
tinuously for U > [/ci, indicating a second-order transi- 
tion from the paramagnetic metallic phase to the metallic 
phase with striped magnetic order. At U > Uc2^ it drops 
abruptly to a small value, corresponding to the formation 
of a full gap in the insulating phase (the finite DOS at 
large U is simply caused by the artificial broadening of 
delta functions to plot results). 



7. Stability of the striped phase 

Consider now the stability of the striped spin-order 
state by analyzing two possible magnetically ordered 
states, one with q=(0,7r) and another one with (7r,7r). 
The U dependence of the magnetization and the en- 
ergy difference between these two states are presented in 
Fig. HQI The (tt, tt) order appears at /7 ~ 2.5 eV, higher 
than the Ud for the (0, tt) state. However, with increas- 
ing [/, the (tt, tt) staggered magnetization increases much 
faster than for the striped state. As shown in the inset 
of Fig. [TOl the (tt, tt) state becomes more stable than the 
(0, tt) state for U ^ 2.65 eV. We have analyzed several 
other sets of hopping parameters and ratios J/U and the 
results are qualitatively the same: the spin stripes are 
stable, but in a narrow region of couplings. 

Hence, in our model the stripe-order ground state be- 
comes unstable to another antiferromagnetic state with 
q=(7r,7r) at ^ 2.65 eV through a first-order transi- 
tion. This is not too surprising considering the values 
of the hopping parameters: the NN and NNN hoppings 
are similar in magnitude. Thus, there is a competition 
between different spin tendencies, as it occurs at inter- 
mediate couplings in Heisenberg spin systems with NN 
and NNN terms. However, the neutron scattering exper- 
iments show that the ground state of several undoped 
Fe-pnictides presents only the striped spin order. To un- 
derstand this predominance of one state over the other, 
note that it has been argued that the emergence of the 
striped state is closely related to a structural phase tran- 
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FIG. 9: (Color online) (a) Density of states at various values 
of U showing the development of a pseudogap at the chemical 
potential in the spin-striped (0,7r) state, (b) The evolution of 
N{co), the DOS at the chemical potential, indicates three well 
separated regions. The non-zero values of N((jJ — /j) in the 
insulating phase at U > Uc2 arises from the finite broadening 
of the raw DOS numerical data. 



sition^ii^ from space group PA/nmm to P112/n. The 
lattice distortion at the transition breaks the four-fold 
rotational symmetry and lifts the degeneracy of the xz 
and yz orbitals. Thus, according to recent calculations 
that incorporate the lattice distortion, the system prefers 
the (0, tt) state over the (tt, tt) state^^ Since the effects 
of lattice distortions are not included in our model, then 
the (tt, tt) spin-ordered state apparently strongly com- 
petes with the spin stripes, but this is misleading and 
caused by the absence of lattice energetic considerations. 

This discussion leads us to believe that simply analyz- 
ing the spin-stripe state found in the mean-field approx- 
imation should be sufficient to understand some of the 
electronic properties of the real system. However, it is re- 
markable that even in the competing (tt, tt) ordered state 
there are also two magnetically-ordered phases: a metal- 
lic phase for moderate U values, and an insulating phase 
with a finite gap at larger U values (see Fig. [TT]) . This 
qualitatively agrees with the previous analysis for the 
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FIG. 10: (Color online) Main panel: the evolution of the 
magnetization with U for q=(0,7r) and (tv^tv). Inset: Energy 
difference between the two states of the main panel. 



striped ordered phase. This suggests that the existence 
of a metallic magnetically ordered phase at moderate U is 
an intrinsic property of the multi-orbital Hubbard model, 
treated in the mean-field approximation, that is robust 
varying the interactions, as discussed in more detail in 
the next subsection. 



with several active orhitals 

We have applied the mean-field technique to several 
other multi-orbital models such as a five-orbital model^ 
another four-orbital model^ and an effective three- 
orbital model including the xz, yz^ and xy orbitalsi^ 
For all these multi-band models a Coulombic interaction 
term similar to Eq. (p!3|) has been used. A robust conclu- 
sion of our mean-field analysis is that a transition from a 
paramagnetic metal to a metallic striped spin-order state 
is found in all the models considered here. 

Let us discuss the results for the five-orbital model*^ 
This model does not give precisely the correct Fermi sur- 
face topology in the undoped case, since it contains a 
pocket at (7r,7r). This problem can be fixed by slightly 
modifying the electronic concentration to, e.g., n = 6.2. 
However, in our study we will consider n = 6.0 for consis- 
tency with the rest of the analysis. Within the mean- field 
approximation, the magnetization at n = 6.0 is shown in 
Fig. [121 In this case there is a broad region of metallicity, 
and the second critical Hubbard coupling (also indicated) 
does not involve a discontinuity. The Fermi surface in the 
inset shows pockets at (0, 0), (tt, 0), (0, tt), and (tt, tt). We 
conclude that a variety of multi-orbital models show sim- 
ilar features as the four-orbital model analyzed before, 
particularly regarding an intermediate metallic magnetic 
phase. 
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FIG. 12: (Color online) Main panel: the evolution of the mag- 
netization for the state with q=(0, tt) in a five-orbital model, 
at electron filling n — 6.0. Inset: the Fermi surface at [7 = 0. 
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FIG. 11: (Color online) Mean- field band structures of the 
spin-ordered (7r,7r) state in the energy window [-2 eV, 2 eV]. 
The upper panel is obtained at JJ — 2.75 eV where band 
overlap indicates a metallic state, even if the order parameter 
is nonzero. The lower panel is at [/ = 4.50 eV where the gap 
is fully developed. 



8. Mean-field results for other models 



III. RESULTS FOR THE TWO-ORBITAL 
MODEL 

The complexity of four- and five-orbital Hamiltonians 
leads to very large Hilbert spaces even for small clusters 
and, as a consequence, it is not possible to compare the 
mean field results against Exact Diagonalization (ED)M 
or Variational Cluster Approach (VCA)— results. For 
this reason the two-orbital model, which was recently 
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studied with ED and VGA methodsfii is here revisited 
to assess the vahdity of the mean-field approximation. 
Our conclusion is that all three techniques lead to simi- 
lar results for the two-orbital model, lending support to 
our claim for the existence of an intermediate-coupling 
metallic and magnetic state. 

The two-orbital model has been widely analyzed in re- 
cent literature and its derivation and other properties will 
not be repeated here. It is given by 

H20 = + ^int, (22) 
where the kinetic energy Hk in real space is)^i^ 

i,cr 

- t2 Xl(4x,a<^i+*.^.<^ + 4y ji+y,y,a + h.C.) 

i,cr 

+ U Y^{dlx^crdi+x+y,y,a + dl^y^^di^x^y,x,a + h.C.) 
i,cr 

— ^4 ^(<^i>,cr^i+^-y,2/,cr + dl^y^^di-^x-y,x,a + h.C.) 

i,cr 

(23) 

The form of Hk in momentum space was provided in 
Refs. andlssi and is given by Eqs. ([3|)-([5j) of the four- 
orbital model. The Coulomb interaction terms are 

^int = ni,a,T^i,a,i + {U' - J/2) ^ rii^xrii^y 

i,a i 

- 2 J ^Si,^ • Si,^ (24) 

i 



where the notation is the same as for the case of the four- 
orbital model but with a = x^y here denoting the orbitals 
xz and yz. The index /i is a unit vector linking NN sites 
and takes the values x or y. ji is the chemical potential. 
As for the case of four orbitals, the relation U'=U — 2 J 
originating in rotational invariance^^ was used. In the 
last term in Eq. ([24|) . the same rotational invariance also 
establishes that the pair hopping coupling J' must be 
equal to J. As before, the hoppings are determined from 
the orbital integral overlaps within the SK formalism or 
from LDA band dispersion fittings J^i^ Most of the prop- 
erties of this two-orbital model are formally similar to 
those of the four-orbital model and, as a consequence, 
several details of the analysis of the previous section do 
not need to be repeated here. 



A. Mean field approximation 

Following the same procedure used for the four-orbital 
model, here we consider only the mean field values that 
are diagonal with respect to the Fe site, orbital, and spin 
labels such that 

{d\^a,cjdv ,cx' ^cj') = (nc, + ^cos(q.ri)mc,)(^ir(^aa'^^aa'. (25) 

After introducing Eq. ([25|) to decouple the four-fermion 
interactions in Eq. ([24|) . and then transforming into mo- 
mentum space, we obtain: 



+ ■ 



^(uny^ 2(U' - i^)nx \ n\^ya 

k,cr ^ ^ 

J ^{Umx + Jmy){dl^Ji,^^xa + 4+q:.a^kxa) 



kcr 



- i ^{Umy + Jmx){dly^di^^^y^ + dl^^y^di^ya). 

kcr 

(26) 

The four mean- field parameters n^^, n^y, m^^, and ruy 
are determined in the usual way by minimizing the en- 
ergy and by requesting that the system be half-filled. 
The half- filling condition determines that Ux = ^y =0.5 
while the values of rux and rriy are a function of U and 
J. We have observed that varying C/, and for fixed J, 
rria becomes non-zero at a critical value U = Ud where 
a gap opens separating the "valence" and "conduction" 
bands at particular momenta but, overall, there is still 
an overlap in energy of some bands and the chemical 
potential crosses both of them. Thus, the system has 
developed magnetic order but it is still a metal. How- 
ever, the bands no longer overlap when U > Uc2 and, 
thus, a metal insulator transition occurs. These results 
are similar to those obtained using more orbitals in the 
model. 

In Fig. [131 the stripe magnetization m=mx + rriy vs. 
U is shown for J=0, /7/4, and U/8. Panel (a) contains 
the results for the LDA fitted hoppings^ while panel 
(b) is for the SK hopping^ with pda = -0.2. The 
critical coupling Ud is the value of U where m becomes 
different from zero. There is a second critical coupling 
Uc2 that is obtained by monitoring the density of states, 
and it separates a metallic from an insulating regime. 
For completeness, this second coupling is indicated with 
a full circle for each case. In panel (b), the shape of the 
curves is the same for the three values of J investigated, 
and the actual value of Ud mildly depends on J. 

The mean-field density of states (DOS) N{u) is pre- 
sented in Fig. [TUfor some values of U and J=U/A^ and 
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FIG. 13: (Color online) Two-orbital model mean-field calcu- 
lated spin-stripe magnetization m vs. U. Panel (a) is for 
the LDA fitted hoppings,— while panel (b) is for the SK 
hoppings.— The solid red line is for J=0. With a blue dashed 
line the results at J =17/4 are indicated. The dotted black 
fine denotes results for J=U/8. The large dot in each curve 
indicates the value of m at Uc2 where the metal to insulator 
transition occurs (see text). Note the absence of a first-order 
transition at Uc2 as in the four-orbital model. 




FIG. 14: (Color online) Two-orbital model mean-field calcu- 
lated density of states (a) for the LDA fitted hoppings;"^^ (b) 
for the SK hoppings. Panel (a): the solid red line is for 
?7=2.5, blue dashed line for U = 5.0, and dotted black line 
for U = 8.0. Panel (b): red is for [7=0.5, blue is for [7=0.8, 
and black is for [7=2.0. J=[7/4 is used in both panels. 



the two sets of hoppings considered here. The solid curve 
in both panels shows results for [7 < Ud- In this case 
the system is metallic. The dashed curve displays the 
DOS for a value of U in between the two critical points. 
Although the DOS varies continuously as [7 increases, 
it is clear that this regime is qualitatively different: A 
deep pseudogap has developed at the chemical poten- 
tial. The system is still metallic in this regime, although 
likely with "bad metal" characteristics. Finally, the dot- 
ted curve shows the DOS for [7 > Uc2' In this case, there 
is a gap at the chemical potential and the system has 
become insulating. Interestingly, the transition at Uc2 is 
not first order for the two-orbital model, in contrast to 
the case with four orbitals. 
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FIG. 15: (Color online) Two-orbital model mean-field band 
structure along high symmetry directions in the extended 
FEZ. The panels on the left column show results for the LDA 
fitted hoppings,-^^ while results on the right are for the SK 
hoppings.i^ (a) [7=2.5, (b) [7=5.0, (c) [7=8.0, (d) [7=0.5, (e) 
[7=0.8, (f) [7=2.0. In all cases, J=[7/4 and the magnetic 
order wavevector is (7r,0). 



Mean field results for the bands and Fermi 
surface 



The mean-field band structure obtained by solving the 
mean-field self-consistent equations is shown along high- 
symmetry directions in the Brillouin zone in Fig. [151 The 
panels on the left (right) column correspond to LDA fit- 
ted (SK) hoppings. The top row shows the band dis- 
persion for [7 < Uci. The corresponding Fermi surfaced 
(FS) in the extended and reduced Brillouin zones (BZ) 
are shown in Figs. [16] (a) and (c), and they agree with 
previous discussions for the two-orbital model. 

The second row of panels in Fig. [15] shows the band 
dispersion in the interesting regime Ud < U < [7^2- It 
can be observed that gaps have opened along, e.g., the 
direction (0, 0)-(0, tt) but not along other high-symmetry 
directions. The partial gaps remove portions of the orig- 
inal FS and produce the 74(k, cj) disconnected features 
(arcs) at T shown in panel (b) of Fig. [16] Note that 
in Fig. [7] (b), the results for the four-orbital model also 
contained arcs, but they appeared in a more symmetric 
manner, namely with four arcs surrounding the T point. 
Two of those arcs were located in the inner hole pocket 
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FIG. 16: (Color online) Mean-field photoemission Fermi sur- 
face for the LDA fitted hoppings^^ in the unfolded FBZ for 
(a) [/ = 1, and (b) [/ = 3, and in the folded FBZ for (c) 
[7=1, and (d) [7 = 3. The ratio J — U /A was used, and the 
magnetic order wavevector is (7r,0). These results were ob- 
tained via A(k,a;) using the same symmetrization procedure 
and energy window centered at the Fermi energy as in Fig. [71 



and two in the outer hole pocket, while in Fig. [16] (b) 
there are only two arcs around the F point. However, as 
often remarked in the context of the two-orbital model, 
to compare with either experiments or with results of 
models with more orbitals, it is crucial to fold the results 
since by this mechanism the hole pocket at M becomes 
the outer hole pocket at F. Following this procedure, in 
Fig. [16] (d) the folded results are shown and now there 
are four arcs around the F point in good agreement with 
the results of the four-orbital model. We conclude that in 
the interesting intermediate coupling regime, both mod- 
els give similar results upon folding of the extended Bril- 
louin zone. 

The last row of panels in Fig. [15] shows the band dis- 
persion for U > Uc2- Now the upper and lower bands no 
longer overlap. The gap is complete and there is no FS. 
The system has become an insulator, as it can be seen in 
the DOS shown in Fig. [H 

The intensity of the features determining the Fermi 
surface should be calculated using the spectral function 
74(k, cj). Figure [15] only shows the eigenvalues of the 
mean-field study, without incorporating the photoemis- 
sion intensity of each state. It is only when the strength 
of the coupling U becomes very large that the spectral 
weights for all the bands will be equal. Otherwise, some 
of the bands will produce strong FS while others will pro- 
duce only weak magnetically-induced "shadow" features 
that are hard to observe, as already shown in Fig. [16] 
To better visualize the bands induced by magnetic or- 
der, in Fig. [17] the mean-field spectral function A(k, uo) 
is presented along high-symmetry directions in the BZ 
for the two-orbital model with the SK parameters For 
U < Uci, panel (a), the spectral weight resembles the 
non-interacting band structure, i.e. there is negligible 
spectral weight in the magnetically-induced bands. For 
Uci < U < Uc2^ panel (b), the bands become distorted 
and bands of magnetic origin develop particularly at the 
locations in which a gap opens. There are other bands 
still crossing the Fermi energy, thus the system is metal- 
lic. Finally for U > Uc2-, panel (c), the gap is complete 
and the magnetic bands are well developed, i.e., four 
peaks can be observed in A(k, uo) for almost all values 
of k. Figure [18] shows the results for the LDA fitted 
hoppings,^^ at J=/7/4: here a similar qualitative discus- 
sion applies. Other values of J such as J=0 and J = U/8 
(not shown) were also considered, and the results are 
qualitatively the same. 



C. Exact Diagonalization results 

To analyze the qualitative reliability of the mean-field 
results, we have performed ED calculations on finite clus- 
ters. In our previous effort^ it was discussed that due 
to the rapid growth of the Hilbert space with the num- 
ber of sites the largest cluster where the two-orbital 
model can be exactly diagonalized has only A^=8 sites. 
It is a tilted ^/Sx^/S cluster, and when periodic bound- 
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FIG. 17: (Color online) Two-orbital model mean-field spectral 
function along high symmetry directions in the extended FBZ, 
using SK hoppings/^ (a) [7=0.5, (b) [7=0.8, (c) [7=2.0. The 
Hund coupling is fixed to J =11/ A and the magnetic order 
wavevector is (7r,0). 
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FIG. 18: (Color online) Two-orbital model mean-field spec- 
tral function along high symmetry directions in the extended 
FBZ for the LDA fitted hoppings.— The couplings used are 
(a) [7=2.5, (b) [7=5, (c) [7=8. In all cases J=[7/4, and the 
magnetic order wavevector is (7r,0). 



ary conditions are implemented the available values of 
the momenta are k = (0, 0), (±7r/2, ±7r/2), (0, tt), (tt, 0), 
and (tt, tt). This limited set of momenta is not well suited 
to analyze band dispersions in the BZ and the use of 
"twisted" boundary conditions (see below) for the 8-site 
cluster would require too high a computational effort. 
However, via the study of spin correlations it has been 
observed that the spin-striped magnetic order character- 
istic of this model is also apparent in the even smaller 
2x2 cluster. For this very small system, the limited 
number of available momenta can be enlarged by imple- 
menting "twisted" boundary conditions (TBC), namely 
requesting that d{Ni + l)=e^^d(l) where Ni is the num- 
ber of sites along the i=x or y direction in the square 
cluster, and (j) is an arbitrary phase. With these TBC 
the values of momenta allowed are now ki = ^^^j^^ with 
Hi ranging from to Ni — 1. Thus, we can calculate 
the spectral functions for a variety of values of k using 



this TBC approach applied to the 2x2 cluster. While 
the very small size is still a serious limitation, note that 
there are simply no other procedures available to contrast 
the mean-field results against exact results at intermedi- 
ate couplings. Our goal using this limited size cluster is 
merely to analyze if mean-field conclusions stand against 
exact results. 

In Fig. [191 the spectral function 74(k, uj) is presented 
along the main diagonal of the extended BZ for differ- 
ent values of U and with J=[7/4. These data have to 
be compared with the mean-field prediction shown in 
Fig. [iHl We present the results for [7=0 for compar- 
ison and to demonstrate that the correct dispersion is 
obtained in spite of the fact that the cluster is so small. 
The finite values of [7 have been chosen to be in the 
magnetic-metallic region ([7=2.5 and 5.0) and in the in- 
sulating region ([7 =8) according to the mean-field results. 
The main point of this figure is to report the develop- 
ment of bands induced by magnetic order with increasing 
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FIG. 19: (Color online) Two-orbital model spectral function 
along the (0, 0) to (tt, tt) direction in the extended FBZ for the 
LDA fitted hoppings.^^ (a) is for [7=0.0, (b) for [7=2.5, (c) 
for [7=5.0, and (d) for [7=8.0. The Hund coupling is J=[7/4. 
The method is ED and the lattice is 2x2 with TBC. The 
arrows indicate the magnetic bands discussed in the text. 



[7. A representative momentum for these magnetically- 
induced bands is highlighted with an arrow in the figure. 
With increasing [7, the magnetic bands smoothly develop 
(concomitant with the development of spin stripe order 
at short distances, as shown later) and at least along the 
main-diagonal direction in the extended FBZ that should 
occur simultaneously with the opening of a gap. Thus, 
our first conclusion is that the extra weak features in 
the one-particle spectral function predicted by mean-field 
due to the magnetic order do appear in the ED results. 
At large [7 there is no doubt that a substantial gap is 
observed, as in the mean-field approach. We also per- 
formed calculations with other values of J such as J = 
and J = [7/8 (not shown), and the conclusions are qual- 
itatively the same as for J=U/A. Qualitatively similar 
conclusions were reached using the SK hoppingsJ^ 

Let us analyze now other directions in momentum 
space. In Fig. [20l the (0, 0) to (tt, 0) and (tt, 0) to (0, tt) 
directions are investigated at [7=0 and 2.5. The results 
indicate negligible changes along these directions by turn- 
ing on [7 : the system appears to remain metallic. How- 
ever, as shown below, the NN spin-stripe order in this 
small cluster is already robust at [7=2.5. Thus, these 
results are compatible with the concept of a state simul- 
taneously metallic and magnetically ordered. Moreover, 
by monitoring the opening of the complete gap we found 
that the metal-insulator transition occurs at a value of 
[7c2 in good agreement with the mean- field predictions. 
The existence of a Ud is a more complicated issue but 
it can be inferred from the development of the magnetic 
bands in A(k, uj) which also occurs in a range of [7 con- 
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FIG. 20: (Color online) Two-orbital model spectral function 
along the (0, 0) to (tt, 0) and (tt, 0) to (0, tt) directions in the 
extended FBZ for the LDA fitted hoppings.^^ (a,b) is for [7=0, 
while (c,d) is for [7=2.5. The Hund coupling is J=[7/4. The 
method is ED and the lattice is 2x2 with TBC. This figure 
shows that the results barely change between the two values 
of [7, suggesting the survival of a metallic state at nonzero [7. 

sistent with mean- field. 

In Fig. [21] we present the magnetic structure factor 
S for k=(0,7r) which is the value of the momentum 
for which it has a maximum (degenerate with (tt, 0) 
for this small cluster); panel (a) shows results for the 
SK hopping^ while panel (b) is for the LDA fitted 
hoppings.-^ To reduce finite-size effects results for the 
TV = 8 cluster are presented, although the results in the 
2x2 cluster are qualitatively similar. For large [7, the 
monotonic increase of S with [7 agrees with the mean- 
field results. At small and intermediate [7, the spin stripe 
correlations are robust and for these small clusters this 
is equivalent to long-range order. But the apparent lack 
of a gap at [7=2.5 along particular directions in momen- 
tum space (discussed before) lead us to believe that the 
ED results are compatible with a metallic and magnetic 
phase at intermediate [7. 



D. VGA results 

In this final subsection, numerical results for the spec- 
tral functions and the density of states obtained with the 
Variational Cluster Approximation (VGA) techniqu o^^i^^ 
are presented. This method embeds the ED solution of 
a small 2x2 cluster into a very large system of a size 
comparable to the 100 x 100 momentum points used in 
the mean field, and thus interpolates between the results 
obtained independently by the ED and mean-field ap- 
proaches. The VGA results discussed here are for the SK 
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FIG. 21: (Color online) The magnetic structure factor at 
k=(0,7r) vs. U calculated by ED of an 8-site cluster for the 
values of J indicated, (a) is for the SK hoppings^^ and (b) 
for the LDA fitted hoppings.^ The insets show the results in 
a more extended range of U. 



parameters. Figure [22] shows the VGA density of states. 
The behavior with increasing U is remarkably similar to 
that obtained in our mean-field calculations presented in 
Fig. [111(b). Metallic, pseudogap, and insulating regimes 
can be clearly observed. 



• u=o.o 




CO-jLl 



FIG. 22: (Golor online) VGA calculated density of states for 
different values of U in the metallic (good metal and pseudo- 
gap) and insulating regimes with the SK hopping parameters 
using pdn/pda^— 0.2 and J^U /A. 

The results for A(k, uo) calculated with VGA are shown 
in Fig. [23l Once again, a remarkable quantitative agree- 
ment with the mean-field results of Fig. [T71 is found. 




FIG. 23: (Golor online) VGA calculated spectral functions 
along high symmetry directions in the BZ for the SK hop- 
ping parameters with pdn /pda=— 0.2, J=U/4, and wavevec- 
tor (7r,0). (a) [7=0.5; (b) [7=0.8; (c) U=2.0. 

IV. CONCLUSIONS 

In this investigation, the mean-field technique was 
applied to multi-orbital Hubbard models for the Fe- 
pnictides. Varying [7, three regions were observed. At 
small coupling, the results are as in the non-interacting 
limit. In the other extreme of very large [7, the ground 
state has a robust gap and the magnetic spin-stripe or- 
der parameter is large. The main result of our effort 
is the presence of an intermediate U coupling regime 
where magnetic spin-stripe order is shown to be compat- 
ible with a metallic ground state due to band overlaps. 
This state has similar characteristics as the parent com- 
pounds of the Fe-pnictide superconductors. Although 
further theoretical work is still needed to firmly estab- 
lish the existence of this interesting intermediate coupling 
regime, for the case of two orbitals our conclusions were 
tested using the ED and VGA methods, and the results 
are compatible with mean-field. 

The analysis of the intermediate U regime allowed us 
to predict the results for angle-resolved photoemission 
experiments. Interesting anisotropics manifested as arcs 
at the Fermi surface. New bands of magnetic origin were 
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also discussed. The two- and four-orbital models lead 
to similar results in this context. Future work will ad- 
dress the optical properties of the intermediate coupling 
regime, and the superconducting state that may arise 
from its doping. 
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